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ABSTRACT 

By considering clans as genuine elementary subprocesses, i.e., in- 
termediate parton sources in the Simplified Parton Shower model, 
a generalized version of this model is defined. It predicts analyti- 
cally clan properties at parton level in agreement with the general 
trends observed experimentally at hadronic level and in Monte 
Carlo simulations both at partonic and hadronic level. In par- 
ticular the model shows a linear rising in rapidity of the aver- 
age number of clans at fixed energy of the initial parton and its 
subsequent bending for rapidity intervals at the border of phase 
space, and approximate energy independence of the average num- 
ber of clans in fixed rapidity intervals. The energy independence 
becomes stricter by properly normalizing the average number of 
clans. 
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Introduction 

Clan concept was introduced in [1] in order to interpret the wide occurrence 
of NB regularity in terms of a two step cascade process. Clans are denned, as well 
known, as group of particles of common ancestor; they are independently produced 
and each clan contains at least one particle. With these assumptions the average 
number of clans is distributed according to a Poisson distribution. Negative Bi- 
nomial (NB) behavior is obtained by requiring in addition that particles inside an 
average clan have a logarithmic distribution which results from the convolution of 
truncated geometric distributions for particles within single clans [2]. 

Clan structure analysis in terms of the average number of clans, N, and of the 
average number of particles per clan, n c , at final hadron level has given striking and 
unexpected results which are still puzzling. In particular it has been shown that the 
density of clans per unit of rapidity in symmetric rapidity intervals, although dif- 
ferent in e + e~ annihilation [3] from deep inelastic scattering[4] and hadron-hadron 
collisions [5], is approximately independent of the cm. energy of the collision; fur- 
thermore the average number of clans at fixed cm. energy grows initially linearly 
with the width of the rapidity interval, then a characteristic bending is seen at the 
border of phase space. 

It should be pointed out that these regularities are satisfied with a higher 
degree of accuracy at final parton level than at final hadron level as it was seen 
by applying clan structure analysis to the final parton Multiplicity Distributions 
(MD's) in symmetric rapidity intervals reconstructed via Generalized Local Parton 
Hadron Duality (GLPHD) hadronization prescription [6] from the corresponding 
MD's at final hadron level [7]. The conclusion was that a complex phenomenon 
like Multiparticle Production might very well be at parton level more elementary 
than at final hadron level. 

The natural question at this point was the following: is simplicity maintained 
at an even more elementary level of investigation, i.e., at single jet level? Single 
quark- and gluon-jets can be isolated as well known in experiments and Monte 
Carlo simulations using a convenient jet-finding algorithm. The use of these algo- 
rithms, although still under discussion, should be considered at present the only 
way which we have to localize single jet properties from many jet events samples. 
In this connection, it was interesting to find in Monte Carlo simulations that the 
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mentioned regularities for clans are valid with a high degree of precision for final 
hadrons in single jets originated by an initial quark and gluon[8]. The accuracy is 
again improved by going via GLPHD to final parton level. 

In parallel, the theoretical study of single parton shower was progressing in the 
literature[9,10,ll,12]. A contribution along this line of thought is the Simplified 
Parton Shower (SPS) model[13], which describes final parton MD's in a single 
jet originated by an initial parton assuming that partons are controlled by the 
dominant interaction among gluons (gluon self-interaction) in a correct kinematical 
framework. The model was solved numerically and results confirm the relevance 
of clan structure analysis for discussing final parton MD's properties. 

New perspectives on the theoretical basis of clan concept were opened later 
on by the discovery of its link with void probability in rapidity space, i.e., the 
probability to detect no particles in a given domain of phase space, and the related 
hierarchical structure of the corresponding correlation functions [14]: the relevance 
of clan concept emerged fully within the class of Infinitely Divisible Distributions 
which reduce to NB MD in a particular case, i.e., the clan concept was found to 
be more general than NB regularity. 

All these undeniable successes notwithstanding, our search failed in explaining 
the independence of the average number of clans on cm. energy in a fixed rapidity 
interval, its linear dependence on rapidity at fixed cm. energy and its subsequent 
bending for larger rapidity intervals. To our knowledge, it should be added, these 
facts are neither explained nor predicted in single parton showers (jets) by any of 
the existing models in the field. 

In order to answer these challenging questions at parton level, we decided 
to implement clan concept in the SPS model[13]. This generalized version of 
SPS, which we call Generalized Simplified Parton Shower (GSPS) model goes 
as follows. We assume that the number of clans in each event coincides with the 
number of steps of the Markov cascade originated by an initial parton, i.e., at each 
step in the evolution of a single jet independent intermediate shower sources, the 
clans, are emitted with degrading virtualities. Consequently, the number of clan 
ancestors is distributed according to a shifted Poisson distribution. By assuming 
local fluctuations in virtuality of the initial parton at the origin of the cascade and 
of clan ancestors, we calculate analytically the energy and rapidity dependence of 
the average number of clans. It is shown that clan properties at partonic single 
jet level obtained from experimental data and from Monte Carlo simulations via 
GLPHD are correctly reproduced by our model. A new striking regularity for the 
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average number of clans is also discovered. 

In section II we review the structure of the SPS model for single jets; in 
section III we present its generalized version and calculate clan MD's both in full 
phase space and in symmetric rapidity intervals. Conclusions are drawn at the 
end. 

II. Summary and critical reading of SPS model for single jets. 

In this section we examine the main features of SPS model which are relevant 
for our scope and which we developed in [13]. Our aim was to obtain the virtu- 
ally evolution and the rapidity structure of a single parton shower by assuming 
essentials of QCD and the validity at each step of the cascading process of the 
energy-momentum conservation law. Inspired by the criterium of maximum sim- 
plicity, essential features of QCD which we decided to incorporate in the model 
are the dominance of gluon self-interaction as described by a kernel of the form 
suggested in [15] and a Sudakov form factor term in order to normalize the ele- 
mentary virtuality splitting function. Virtuality evolution in a single jet (parton 
shower) was then described as follows. 

We consider an initial parton of maximum allowed virtuality W which splits 
at virtuality Q into two partons of virtuality Qo and Q\. We require Q > Qo + Qi 
and Qo, Qi > 1 GeV. We define the probability for a parton of virtuality W to split 
at Q, p(Q\W), which is normalized by a Sudakov form factor. The probability for 
an ancestor parton of maximum allowed virtuality W to generate n final partons, 
P n {W), and the probability for a parton which splits at virtuality Q to generate n 
final partons, R n (Q), with R n (Q) = S n i for Q < 2 GeV, lead to the corresponding 
generating functions: 

CO 

f(z,W) = J2Pn(W)z n - 1 (1) 
n=l 

CO 

g(z,Q) = J2Rn(Q)z n - 2 (2) 

n=2 

The two generating functions are linked by 

f(z,W) = J p(Q\W)dQ + J p(Q\W)zg(z,Q)dQ (3) 

The joint probability density V(QoQi\Q) for a parton of virtuality Q to split into 
two partons of virtuality Qo an d Qi is defined by 

V(QoQi\Q) = p(Qo\Q)p(Qi\Q)K(Q) (4) 
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where K(Q) is a normalization factor. 

The dynamical content of the model in virtuality is contained in the following 
equation 



n-l 



" poo />oo 

Rn{Q)=Y J dQ dQ 1 V(Q Q 1 \Q)R n - n <(Qo)R n >(Qi)e(Q-Qo-Qi) 

n' = l Jl Jl 

(5) 

which gives the probability for a parton which splits at virtuality Q to generate n 
final partons in terms of the joint probability density for a parton of virtuality Q 
to split into two partons of virtuality Qo an d Qi- 

Eq. (5) can be reformulated for the corresponding generating function by dividing 
the domain of integration in three subdomains; one obtains 

dQ J dQ 1 V(Q Qi\Q)e(Q-Qo-Qi) + 

/I poo 
dQo J dQ 1 V(Q Qi\Q)zg(z,Q 1 )e(Q -Q -Qi)+ (6) 

pOO pOO 

/ dQo / dQ{P{QoQi\Q)zg{z, Q )zg(z, Q 1 )9(Q -Q - Q x ) 

J2 J2 

The above three subdomains correspond to the possible different situations in 
which the two generated partons can be found i.e., no one or only one or both 
split. This general scheme is valid for any splitting function p(Q\W). In case 
p(Q\W) is factorizable in terms of its variables Q and W, i.e., 



p(Q\W)= Po (Q)C(W), C(W) = 

eq. (3) simplifies into the differential equation 

df(z,W) 



w 

dQ'po(Q') 



(7) 



m] p(W\W)[zg(z,W)-f(z,W)} (8) 

For numerical simulations we choose: 

^niwwn A (\ogQ) A -\ ^ /logQ\ A 
p(Q\W)dQ = — {lQgW)A dQ = d j (9) 

A being the only free parameter of the model. This form of p(Q\W) was motivated 
in [13] by our request of simplicity in the structure of the model. Notice that eq. 
(9) corresponds to the virtuality dependence of the standard QCD kernel. 



4 



For describing the rapidity structure of the model, we proposed in [13] to use 
the singular part of the QCD kernel controlling gluon branching: 

p(yo\QoQiQy) oc P(z )dz oc — - (10) 

2o(l - z o) 

Here z$ is the energy fraction carried out by the produced parton in the infinite 
momentum frame. 

The limits of variation of zq are fixed by the exact kinematical relations 



where B = ^[1 4- (^=f) 2 — (^q~) 2 ] is the scaled parton energy in the center of mass 



system and y B 2 — ( g ) 2 its maximum scaled transverse momentum. In this way 
the scaled transverse momentum of the parton with energy fraction zq 



^v B2 -(t) 2 - (2 °- B)2 (12) 



vo = v + h°g ~ ' H" ~' (13) 



and its rapidity 

ii, * = ii 4- 

are uniquely determined. Rapidity of the second parton of virtuality Qi is obtained 
by energy-momentum conservation: 



B + 


\zo 


-B\ 


B- 


\zo 


-B\ 



Vi = V + tanh 1 



B 

tanh \y - y\ 



B + l 



(14) 



Notice that only the first step has to be treated differently in rapidity because 
it corresponds to the degrading from the maximum allowed virtuality W to the 
virtuality of the first splitting Q. In this case the rapidity of the ancestor is fixed 
by conservation laws and is given by 



^tanh" 1 t/l"(f) (15) 

The kinematical structure of the model is summarized in Figure 1. 

The virtuality evolution equations in full phase space have been analytically 
solved in two cases: a Poissonian distribution has been obtained when only one 
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produced parton can split and a geometric distribution is found when we neglect 
the virtuality conservation law in the production process. 

We found that the numerical solution of the master equation for the MD in a 
single jet is well reproduced by a NB distribution in full phase space, in symmetric 
rapidity intervals and in p T intervals. It is also found that the average number of 
clans approximately scales with jet energy in a fixed rapidity interval. 

At this level of investigation, it should be clear that still open problems are 
the lack of the analytical solution of equation (6) and, in more general terms, the 
lack of a complete analytical study of the parton evolution process in rapidity and 
p variables. 

T 

III. The Generalized Simplified Parton Shower Model (GSPS) 

In order to solve part of the problems indicated at the end of the previous 
section, we propose to incorporate in the SPS model the clan concept; we call 
this version of the model Generalized Simplified Parton Shower (GSPS) model. 
Accordingly, we decide now to pay attention for each event to the ancestor which, 
splitting n times, gives rise to n subprocesses (one at each splitting) (see Figure 
2) and we identify them with clans. Therefore in this model for a single event the 
concept of clan at parton level is no more as it was in [1], i.e., only a statistical 
one. In the present picture the clans are independent active parton sources and 
their number in each event coincides with the number of splittings of the ancestor, 
i.e., with the number of steps in the cascade. 

Notice that each clan generation is independent of previous history (it has 
no memory); thus the process is markoffian. Furthermore each generation process 
depends on the evolution variable only and is independent of the other variables 
of the process like the number of clans already present and their virtualities. In 
the original version of SPS the splitting function of the first step, p{Qo\Q), was 
different from the splitting function of all the other steps, [3(Qo\Q), obtained by 
integrating the joint probability function V(QoQi\Q) 

P(Qo\Q) = £ Q ° dQ^iQoQilQ) (16) 

In order to generalize the SPS model as it stands we assume that virtuality con- 
servation law is locally violated (although conserved globally) according to 

l<Qo<Q, 1 < Qi < Q (17) 
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The upper limit of integration of (16) becomes Q, the normalization factor K(Q) 
reduces to 1 and the process becomes homogeneous in the evolution variable since 

rQ 

P{Qo\Q) = p(Qo\Q) J p(Qi\Q)dQi=p(Qo\Q) (18) 

where p(Qo\Q) is factorized according to eq. (7). The approximation described 
by eq. (18), therefore, can be interpreted as the effect of local fluctuations in 
virtuality occurring at each clan emission. 

This violation of the virtuality conservation law spoils of course the validity 
of the energy-momentum conservation law, which, in the SPS model, uniquely de- 
termines the rapidity of a produced parton, given its virtuality and the virtuality 
and rapidity of its germane parton (see eq. (14)). In the GSPS model the two 
produced partons at each splitting are independent both in virtuality and in ra- 
pidity; however, the rapidity of each parton is bounded by the extension of phase 
space fixed by its virtuality and the virtuality of the parent parton: 

\yi-y\< log— (19) 

In conclusion, by weakening locally conservation laws, we decouple the produc- 
tion process of partons at each splitting. Consequently, the GSPS model allows 
to follow just a branch of the splitting, since each splitting can be seen here as 
the product of two independent parton emissions. This consideration will be par- 
ticularly useful in discussing the structure in rapidity of the model; in fact, it is 
implied that the Altarelli-Parisi kernel given in eq. (10) should be identified with 

p(yo\QoQy)dy p(yi\QiQy)dy 1 oc (20) 

In the following we propose to study clan formation in a parton shower first 
in virtuality variable and, next, in virtuality and rapidity variables. The p T de- 
pendence of clan properties in the GSPS model could be treated with the proper 
modifications in an analogous way. 

In order to stress the fact that from now on we will pay attention to clan 
production, we use lowercase letters (r,p) when we refer to the probability of 
producing N clans, whereas in section II we used uppercase letters (R,P) for the 
probability of producing n partons. In addition, in order to be extremely clear, we 
proceed to calculate first the probability for emitting N clans without including 
the first step, rjv, and later the same probability by including it, pn- 
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III. a. Clans formation in a parton shower: Full Phase Space 



Let us consider the probability r;sr(Qo\Q) that an ancestor parton, splitting 
at virtuality Q > 2 GeV, gets virtuality Qq after emitting N — 1 clans. 

In the GSPS model this probability can be expressed in terms of the elemen- 
tary splitting function, p(Qo\Q), (see Figure 3, solid line diagrams): 

ro(Qo\Q) = 0, ri(Qo|Q) = 0, r 2 (Q \Q) = p(Qo\Q), 
"Q 



r 3 (Qo\Q)= I dQ lP (Q \Qi)p(Qi\Q), 

J max{2,Q } 



rQ pQn-2 

rN(Qo\Q) = / dQ N - 2 / dQ N - 3 ... 

«/max{2,Qo} «/max{2,Qo} 

-Q2 

dQ lP (Q \Qi)p(Qi\Q2) ■ ■ -p(Qn-2\Q) 

'max{2,Q } 



(21) 



J m 



Notice that, since a parton of virtuality less than 2 GeV does not split, the prob- 
ability r/v(Qo|Q) has a different form when the virtuality Qo is larger or smaller 
than 2 GeV. By taking into account the factorization of p(Qo\Q) given in eq. (7), 
eq. (21) can be rewritten as: 



rQ 

(Qo\Q) = P(Qo\Q) / p(QN-2\Q N -2)dQ N -2 

(ZZJ 



/max{2,Q } 

pQn-2 rQ2 

/ p(Q N - 3 \Q N -3)dQ N - S ... p(Qi\Qi)dQi 

</max{2,Qo} </max{2,Qo} 

The multiple integral of the product of identical factors over the ordered domain 
max{2,<5o} < Qi < Q2 < • • • < Qn-2 < Q can be solved by extending the 
integration domain to a (N-2)-dimensional hypercube of length Q — max{2,<5o}- 
Then, by using the symmetry of the integrand, one gets 

(n]n) . P(Qo\Q) [Xi % X ^r 2 Qo>2GeV 
rN(Qo\Q) = { _ , _ [A(Q)] iv-2 ; „ (23) 



p(Qo\Q) 7^ 2) , Q <2GeV 



with 



X(u) = p{u'\u')du' 



The probability tn{Q) that a parton, splitting at virtuality Q, generates a shower 
with N clans is given by integrating Q over the range allowed for a final parton: 

r N (Q) = J dQ r N (Q \Q) (24) 

Since only the second of eq. (23) contributes, one obtains a shifted-Poisson distri- 
bution 

r O (Q) = r 1 (Q)=0 

M0) = N> _ 2 ™ 

The average number of clans generated by an initial parton splitting at Q follows: 

00 00 \\(n\i N - 2 

N r (Q) = £ Nr N (Q) = Ne ~ HQ) (N -2V = HQ)+2 (26) 

N=0 N=2 \ >' 

Now we can take into account the first step and study the multiplicity distribution 
for clans produced from an ancestor of maximum allowed virtuality W (see Figure 
4, solid line diagrams). The probability that an ancestor parton of maximum 
virtuality W gets virtuality Qo after emitting N — 1 clans, Pn(Qo\W), is given by: 

Po(Qo\W) = 0, pi(Q \W) = p(Qo\W), 
p 2 (Q \W)= / dQp(Q \Q)p(Q\W), 

J max{2,Q } 

: (27) 

rW rQ 

Pn(Qo\W) = dQ dQ N - 2 ... 

Jmax{2,Qo} </max{2,Qo} 

<J<9ip(<9oIQi)p(Qi|Q 2 ) ■ ■ .p(Q N -2\Q)p(Q\W) 

'max{2,Q } 



J m. 



It should be pointed out that eq. (27) for N > 2 can be obtained also by rewriting 
Pn(Qo\W) in terms of r N (Q Q \Q): 

Pn(Qo\W) = / dQp(Q\W)r N (Q \Q) (28) 

J max{2,Q } 

This relation for clans corresponds to eq. (3) for partons. 
Explicitly, 

fnmn , P^l^ ^-ff"" 1 Qo > 2 GeV 
Pn(Qo\W)=< ,w m U-i (29) 

piQolW) ^^ Qo < 2 GeV 

9 



Then, by integrating over Qq in the allowed virtuality range, one gets the multi- 
plicity distribution p N (W): 



Po(W) = 

p „ W ,/«M W ^ e ->'»« (30) 

Being included in eq. (30) the graph corresponding to the case in which the 
ancestor does not split (one-parton shower), N starts from 1 and not from 2 as in 
eq. (26). 

Eq. (30) is a shifted Poissonian distribution with average number of clans 
given by: 

N(W) = N Pn{W) = Ne~ x ^ M_JJ = A(W) + 1 (31) 

N=0 N=l ^ >' 

which is fully consistent with eq. (26) via eq. (28). We stress that this result has 
been obtained a priori in the present generalized version of the model, differently 
from what has been done previously in [1] where the independent production of 
clans was an "ansatz" introduced a posteriori in order to explain the occurrence 
of NB regularity. 

Eqs. (30) and (31) answer to our main questions in full phase space. In 
order to introduce the subsequent discussion on rapidity dependence, we propose 
to study first the clan density in virtuality, i.e., the number of clans emitted with 
virtuality Qq from an ancestor parton splitting at virtuality Q. 

We start by considering the probability that the (N — l) th clan is emitted 
with virtuality Qq. Now, it should be remarked that in the GSPS model the 
elementary splitting function for a parent parton of virtuality Q to produce a 
clan of virtuality Qi, p(Qi\Q), has the same functional form of the probability to 
produce the associated parton of virtuality Qq, p(Qo\Q). Therefore, the role played 
by the ancestor and the clan at the last step of the cascade can be interchanged; 
the probability that an ancestor, splitting at Q, gets virtuality Qq after emitting 
N — 1 clans, r/v(Qo|Q)> (see Figure 3, solid line diagrams and eqs. (21,23)) turns 
out to be equal to the probability to produce the (N — l) th clan with virtuality 
Qq (see Figure 3, dashed line diagrams). Finally, by summing over N , the average 
number of clans produced with virtuality Qq, r^(QQ\Q), is obtained; the explicit 
expression of this probability is quite subtle, since one has to carefully take into 
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account the role of the ancestor. For Q > 2 GeV, the ancestor will generate other 
partons and will not contribute to the final clan density: 

OO 

r s (Q |Q) = r ^(Qo\Q) Qo > 2 GeV (32) 

N=2 

For Q < 2 GeV, the ancestor is considered as one parton clan and will contribute 
to r^(Qo\Q) in addition to the term given in eq. (32); thus one has 

oo 

r E (Q |Q) = r N(Qo\Q) + r anc (Qo|Q) Qo < 2 GeV (33) 

N=2 

where r anc (Qo|Q) is the probability to obtain an ancestor of virtuality Qo < 2 
GeV in any number of steps and is given by 

oo 

r an c(Qo|Q)= ^2r N {Qo\Q) (34) 

N=2 

Notice that to eq. (34) will contribute just the second of the two eqs. (23). 
Accordingly 

r a „c(Qo|Q) = £ P(Qo\Q) W=^ = P(Qo\Q)e X(Q) (35) 

N=2 \ >' 

and the full expression for the clan density in virtuality, r E (Qo|Q), turns out to 
be 

r (On\0) ~ I P(Qo\Q)e HQ) - X(Qo) Qo > 2 GeV . , 

r E (Qo|Q) - | 2p{QolQ)em) Qq < 2 GeV (36) 

The consistency of this result can easily be checked: by integrating r s (Qo|Q) over 
Qo i n the full domain one should obtain indeed for the average number of clans, 
iV r (Q), the same expression which was calculated in eq. (26). In fact 

rQ z-2 rQ 

/ dQ r E (Q |Q)= / dQ 2p(Qo\Q)e x{Q) + / dQ oP (Qo\Q)e x{Q) - x{Qo) = 

Jl Jl J2 

= 2 + A(Q) = N r (Q) 

(37) 

The same procedure can be used for studying the density of clan produced by an 
ancestor of maximum allowed virtuality W. We start by considering the probabil- 
ity to produce the (N — l) th clan at virtuality Qo from the ancestor of maximum 
allowed virtuality W (see Figure 4, dashed line diagrams). In the model this prob- 
ability turns out to be equal to the probability that the ancestor of maximum 
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allowed virtuality W gets virtuality Q after emitting N — 1 clans, i.e., to the 
probability Pn(Qo\W) given by eq. (29) (in Figure 4, solid line diagrams). 

Notice that the showers with one parton only do not spoil our reasoning; how- 
ever, they will contribute to the average number of clans produced with virtuality 
<2o, p^iQolW), for Q < 2 GeV. It follows 

v (Qnm _ { Zn=2Pn(Qo\W) Qo > 2 GeV 

where p a nc(Qo\W) is given by definition by 

P,nc(Qo\W) = PN(Qo\W) = P(Qo\W) [ { " 

N=l N=l ^ >' (39) 

= p(Q \W)e x ^ 

The sum starts of course from N = 1 since, as we discussed previously, we must 
include the contribution of one-parton showers. 
Accordingly, one has: 

V (Qa\W) ~ l PiQolW) [e HW) - X(Qa) - 1] Qo > 2 GeV 

The consistency of our calculation can be checked by integrating eq. (40) over Qo 
from 1 to W; eq. (31) follows. 

Eq. (31) solves our first problem to determine analytically the average number 
of clans generated by an initial parton of maximum allowed virtuality W in full 
phase space. It should be added that eq. (36) and eq. (40) are simply connected 
by the relation 

P (Qo\W) = I fZ d ^ W K^\Q) Qo > 2 GeV 

E \C dQp(Q\W) rj: (Q \Q)+p(Qo\W) Q <2 GeV 

which is obtained by introducing eqs. (28), (32) and (39) into eq. (38). 

In the following subsection, eq. (41) will be extended to include rapidity 
variable in order to calculate the average number of clans in symmetric rapidity 
intervals. 

IILb. Clans formation in a parton shower: Symmetric rapidity intervals 

In this subsection we extend our study to clan distributions in symmetric 
rapidity intervals Ay = [— y c ,y c ]- We relate the probability to have N' clans in 
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the symmetric rapidity interval Ay, PN'(y c , W), to the corresponding probability 
defined in full phase space (fps), pjv(W), by using the following relation 

oo 

p N ,(y c ,W) = n(N>,y c \N,fps)p N (W) (42) 

N=N' 

where H(N' , y c \N, fps) is the conditional probability to have N' clans in Ay when 
one has N clans in full phase space; it contains all dynamical information on the 
rapidity structure of the production process. 

Following the discussion in subsection III. a and in particular the fact that 
clans are not correlated in full phase space and the only allowed correlations are 
among particles within the same clan, the conditional probability IL(N' ,y c \N, fps) 
turns out to be a positive binomial distribution 

U(N',y c \N, fps) = (^^'(l -n) N - N ' (43) 

where 7r(y c , W) is the probability to produce a single clan in the rapidity interval 
Ay from the maximum allowed virtuality W. In terms of generating functions 
/clan (2) = J2n=o Pnz n , the above relation can be written as 

/cfa y nW = /c f fa S n(^ + l-^) (44) 

Being clans distributed in full phase space according to a shifted Poisson distri- 
bution (eq. (30)), via eq. (44), the probability generating function for clans in 
symmetric rapidity intervals has the following form 

/ctanW = [7r(y c ,W)z + l-7r(y c ,W)}e x ^> w K*-V (45) 

where A(y c , W) = n(y c , W)X(W). 

Eq. (45), being the sum of two Poissonian distribution (the first a shifted 
one), has a nice physical meaning: the two terms correspond to the probability 
of having the ancestor within or outside the given rapidity interval. Since in full 
phase space 7r(y c = fps, W) = 1, from eq. (45) one obtains the generating function 
for full phase space. In addition, it should be pointed out that the void probability 
in a given rapidity interval, po(y c , W), is different from zero and is given by 

p (y c , W) = [l- AVc W)\ e -*<Ve.w)A(W0 ( 46 ) 

When y c is very small, 7c(y c , W) tends to zero and the unshifted Poissonian dom- 
inates. Thus, it can be stated that in the smallest rapidity intervals the clan 
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multiplicity is to a good approximation Poissonian and the full MD belongs to 
the class of Compound Poisson Distribution[14]. 7r(y c , W) tends to 1 for y c at the 
border of phase space and the exact full phase space shifted-Poisson distribution 
is approached. Finally, when \(y c , W) is sufficiently large, the shifted-Poissonian 
dominates at large N (the tail of the distribution) while the unshifted one dom- 
inates at small N (the head of the distribution). These facts might have some 
consequences in interpreting the anomalies found in NB behavior for small N and 
the deviations from NB behavior in large rapidity intervals, which are controlled 
by the behavior of the distribution at large N. 

From eq. (45) it follows that the average number of clans is given by: 

N(y c ,W) = 7r(y c ,W)N(W) (47) 

We are ready now to calculate the average number of clans in symmetric 
rapidity intervals. 

The probability that the ancestor, splitting at virtuality Q and rapidity y, 
gets virtuality Qq and rapidity yo emitting N — 1 clans, rN(Qoyo\Qy), can be ex- 
pressed within the GSPS model (as for rN(Q \Q)) in terms of elementary splitting 
functions according to the following formula (see Figure 3, solid line diagrams): 

rQ rQ2 

rN(Qoyo\Qy)= dQ N _ 2 ... d<2ip(<2ol<3i) • • -p(Qn-2\Q) 

«/max{2,Qo} «/max{2,Qo} 

/oo />oo 
dy N - 2 ... / dyxYdyo - yi\,Qo,Qi) . . .Y(\y N _ 2 - v\,Qn-2,Q) 
-oo J — oo 

(48) 

where Y(\yi — | , Qi, Qi+i) is the probability that a parent parton with rapidity 
yi + i generates in a single step a parton of rapidity y^, being their virtualities Qi+i 
and Qi respectively. 

Inspired by the criterium of simplicity which we decided to follow from the 
beginning of our work and by the considerations discussed at the beginning of 
Section III, we take for the probability Y(\yi — yi+\ |, Qi, Qi+i) the simplified 
QCD-splitting kernel 

dzi 

Y{\Vi - Vi+i\,Qi, Qi+i)dyi = P(zi)dzi oc — (49) 

Zi 

By using eq. (13) properly rewritten, the kernel becomes 

y{\Vi - Vi+i\,Qi,Qi+i)dyi oc dyi (50) 
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which after the normalization in the kinematically allowed domain \y t — < 
log(Qi + i/(5i) turns out to be: 

Y{\yi - y i+ i\,Qi, Q i+ i)dyi = ^ — -9(\og(Q i+1 /Q % ) - \y t - y i+1 \) (51) 

2 log ' Qh ' 



Notice that in the limit of Qi+i — > Qi the function Y(\yi — Qi, Qi+i) reduces 
to a Dirac 6- function; moreover, by inserting eq. (51) into eq. (48) and by inte- 
grating over yo in all the allowed domain one recovers immediately the full phase 
space result given in eq. (21). 

In this equation, as well as in the case of full phase space, one can identify the 
probability rN(Qoyo\Qy) with the probability that the (N — l) th clan is produced 
with virtuality Qo and rapidity yo from an ancestor which splits at virtuality Q 
and rapidity y (see Figure 3, dashed line diagrams). 

The integrations in the allowed domain in rapidity of eq. (48) are awkward 
and quite difficult. The idea is to approximate the probability density given by 
eq. (51) with a gaussian of the following form: 

YQyi - yi+i\,Qi,Q i+ i) = — L- exp (- \ y * 2 ^ +1 ^ ^ (52) 
of being the width of the distribution given by 

"'M 9 ^)- 10 ^) <53) 

This approximation corresponds to weaken locally energy-momentum conservation 
laws and simulates the onset of local fluctuations on the probability density defined 
by eq. (51). It is remarkable that with the just mentioned simplification the 
integrations on the rapidity domains can easily be performed and lead to the 
following result: 

/OO />oo 
dy N - 2 ... j dy\Y(\yo-yi\,Q Q ,Q{)...Y(\y N -2 -v\,Qn-i,Q) = 
-oo J —oo 

(54) 



1/2 eXP \ 9 9 



\Vi ~ Vi+1\ 



From eq. (53), one sees that the width of the right part of eq. (54) depends only 
on initial and final virtualities, i.e., 

N-2 N-2 . . 

E °i = a E ( lo § - log Qi) = a log ^- (55) 
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It is clear that eq. (48) can be written now in terms of the probability tn(Qq\Q) 
defined in eq. (21) as follows: 



rQ 

rN(Qoyo\Qy) = Y(\y -y\,Q Q) / dQ N 2 

J max -[2. On]- 



max{2,Q } 

dQip(Qo\Qi) . . -p(Qn-2\Q) 

'max{2,Q } 

= Y(\y -y\,QoQ) x r N (Q \Q) 



(56) 



By recalling eq. (23), one has: 

rN(Qoyo\Qy) = 



p(Qo\Q)Y(\y - y\,Q Q Q) [A(Q ^l%, r ~ 2 Qo > 2 GeV 



(AT-2)! 

P (Qo\Q)Y(\y -y\,QoQ) [h $Sf Qo<2 GeV 

(57) 

In summary, the use of gaussian approximation for the step function (step function 
— > gaussian) allows to see that the convolution of a generic number of distributions 
in rapidity domains given by the corresponding step functions is equivalent to a 
single gaussian distribution in rapidity. If we assume that this approximation can 
work the other way around (gaussian — > step function) the initial form of the 
probability density Y{\yi — Qi, Qi+i) as given by eq. (51) is re-obtained 

restoring the correct kinematical domain for the total distribution. 

From eq. (57), it is quite easy to calculate the clan density in virtuality and 
rapidity produced by an ancestor which splits at virtuality Q and rapidity y. Once 
again, in this case the ancestor plays a very peculiar role: for Qo ^ 2 GeV (the 
ancestor does not contribute) one has 



r . 



: (Qoyo\Qy) = Yl r N(Qoyo\Qy) Qo > 2 GeV (58) 

N=2 

and for Q < 2 GeV (the ancestor contributes additively) the probability becomes 

oo 

r^iQoyolQy) = ^ r N(Qoyo\Qy) + r anc (Qoyo\Qy) (59) 

N=2 

where r an c(<2o2/o|<22/) is given by definition by 

oo 

r anc (Qoyo\Qy) = ^2 r N(Qoyo\Qy) (60) 



N=2 
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i.e. 



XQoVolQy) = 2 ^ r N (Q y \Qy) Qo < 2 GeV 



(61) 



N=2 



Finally, we find 



r^QoyolQy) - I MQolQ)Y{ \ yo _ y lQ oQ)e HQ) Qq < 2 G eV (62) 

Clan density in rapidity is calculated by integrating the above density over the 
virtuality Qo i n the full virtuality range: 

rQ 

r^{yo\Qy) = J dQor^iQoyolQy) (63) 
Let us study separately this integral for Q < 2 GeV and Qo > 2 GeV, i.e., 

r^(yo\Qy) = r E ,Qo<2 + ^e,q >2 (64) 

For Qo < 2 GeV one has 

r 2 

rE,Q <2 = 2 J dQ p(Qo\Q)Y(\yo-y\,QoQ)e X{Q) = 



_ 2 A(Q) f 2 dQo AjlogQo^- 1 1 

Ji Qo (logQ) A 21og(£) 



( log ( ^ ) - \yo -y\ 



By defining then logQo/logQ = x, equation (65) can be rewritten as 



rE,Q <2 = A 



(logQ) A_1 /" min ( 1 °g 2 / 1 °g ( 3' 1 -|yo-y|/iogQ) X A ~ X 



(log 2)' 



e&r- 



1 -x 



(65) 



(66) 



Notice that the integral in eq. (66) can be expressed as a series, i.e.: 



ax- = > 

1 -x ^ A + n 

n=0 



(67) 



As an example, let us consider the solution of eq. (66) corresponding to A=2. One 
gets: 



7"£,Qo<2 = 



r S,Q <2 



21ogQ 
[log2] 2 

21ogQ 
[log2] 2 



log 2 
"logQ 



-1 + 



log 1 



log 2 
logQ, 

for \y -y\ < logQ - log 2 

\yo -y\ , (\yo-y\ 



logQ 



-log 



logQ 



(68) 



for |y -y| > logQ - log 2 
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Coming now to the second term of eq. (64) (Q > 2 GeV), one has 



r £,Q >2 _ 



dQ p(Qo\Q)Y(\y - y\,QoQ)e x{Q) - x(Qa) = 



2 



« dQ A(\ogQ ) A -i „ A(Q) _ A(Qnl 1 a ( Q 



2 



o n ^ _AWoJ 777T* lo « 7T ) — l2/o — 2/1 = 

f Q AdQp 1 / /Q\ A 

= / 2 Qoioggo 21og (^)% log UJ- |y °- ?/l J 



(69) 



In terms of x, this equation can be rewritten as 



/-i-lj/o-yl/iogQ ^ 

Notice that the integral diverges for x — > 1, ie., for |j/ — J/| = 0. However, this 
singularity is under control: it is originated by our choice to follow just one branch 
of the splitting process, thus allowing Q — > Q. Although the probability density 
is infinite for \yo — y\ =0, the resulting probability is finite for any finite range of 
yo. Therefore, one has 



= 2l^Q l0g 



logQ- 


yo 


- y\ 




yo - 


-y 





log Q - log 2 



log 2 



\yo - y\ > o 

(71) 

The expression of Eq. (71) corresponding to A = 2 is trivial. 

It should be noticed that Eq. (71) is defined in the kinematically allowed 
region 

0< \y -y\ < logQ-log2 , (72) 

to be compared with the kinematically aloowed region of Eq. (65) (\yo— y\ < logQ). 

In conclusion, by adding eqs. (65) and (71) one obtains r^yolQy), the number 
of clans of rapidity yo generated from an ancestor splitting at virtuality Q and 
rapidity y. The density of clans generated from an ancestor of maximum allowed 
virtuality W, p^(yo\W), can be studied by integrating over Q in the region shown 
in Figure 5 (see eq. (41) where analogous calculations are performed in full phase 
space). The solid line inside the (yo, Q)-domain in Figure 5 corresponds to eq. 
(15). It follows 

rW 

P^yo\W) = J dQp(Q\W)r^yo\Qy)+ Pfmc (y \W) (73) 
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Pa,nc(yo\W) being the probability that the ancestor does not split and has rapidity 
y fixed by eq. (15): 

p anc (yo\W) = J dQ oP (Qo\W)5 (y - tanh" 1 y/l - (Qo/W) 2 ) (74) 

This contribution corresponds to the first diagram of Figure 4 and to the thick solid 
line in the lowest right corner of the (yo, Q) domain of Figure 5. The 6- function 
in eq. (74) can be rewritten in terms of Q$ 

S U - tanh"' v/1 - (Q /W)A = 6 ( Qo - (75) 
V ) \ cosh y J cosh 2 y 

and eq. (74) becomes 

/ | Tin 4 , , [log W - log cosh y ] A_1 
Panc (yo\W) = A tanh y ^ 

for tanh" 1 y/l - (2/W) 2 < y < tanh" 1 y/l - (1/W) 2 

The density of clans in rapidity produced from an ancestor of maximum allowed 
virtuality W can now be calculated by using Eq. (73). Accordingly, the average 
number of clans in symmetric rapidity intervals is given by: 

N(y c ,W)= r dy oPl ,(yo\W) (77) 

J-Vc 

This result is remarkable. It should be pointed out that it depends on parameter 
A only. Parameter A controls - as we have seen in [13] - the number of clans 
produced in full phase space, i.e., the length of the cascade (since the number of 
clans in GSPS model is identified with the number of steps in the shower, more 
steps we have more extended turns out to be the shower). Analytical solutions 
of eq. (77) can be obtained explicitly for all integer positive values of parameter 
A, and numerical solutions can be computed for all positive real values of A. 
We decided to postpone the discussion on the A-dependence of our formula to a 
forthcoming paper and to limit ourselves to indicate the structure of our formulae 
by solving analytically the case A = 2. This limitation does not alter the main 
features of our results; in fact preliminary analysis shows that our results are 
qualitatively not altered by other choices of parameter A. In order to perform the 
integration, we used the approximation 




tanh ' \/1-(t-H ^ log W + log 2 - log Q (78) 
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which slightly modifies the (y , Q) domain (solid line in Figure 5) into the dashed 
lines in the same Figure. We summarize in the Appendix the long and cumbersome 
calculations which are needed in order to integrate eq. (77). 

In Figure 6 the clan density Pj,(yo\W) corresponding to eq. (73) is shown as 
a function of yo variable for maximum allowed virtualities W = 50 GeV, W = 
100 GeV and W = 500 GeV. The contribution of one-parton showers turns out to 
be negligible for this choice of the A value. Notice that the height of the curve is 
decreasing and the width increasing with the energy. Convolution of clan density 
for two parton showers is shown in Figure 7. Notice that the central dip at y ~ 
is slowly removed by increasing the energy of the initial parton. It should be kept 
in mind that the structure of Figures 6 and 7 refers to clan production; found 
different behavior for parton production is not in contradiction with this behavior 
since we have still to include in our scheme parton production within a single clan. 

In Figure 8 the average number of clans N(y c , W) (eq. (77)) is given as a 
function of rapidity width y c for the same W values of Figure 6. Limitations on 
the rapidity intervals are determined by the available phase space corresponding 
to the different initial parton virtualities. 

Accordingly, the GSPS model predicts for the average number of clans 
N(y C: W) at parton level in a single shower (jet): 

a) a rising in rapidity width y c for different initial parton virtualities W very 
close to linear for 1 < y c < yf ps ; the rising is still linear but with a somewhat 
different slope for y c < 1. Characteristic bending occurs finally for rapidity 
width y c < y fps ; 

b) approximate (into 5%) energy independence in a fixed rapidity interval y c for 
W below 100 GeV. For higher virtualities deviations from energy indepen- 
dence become larger; they are into 20% when comparing N(y c , 50 GeV) and 
N(y c , 500 GeV). It should be noticed that the average number of clans slowly 
decreases with virtuality; this behavior has been already observed in Monte 
Carlo simulations for single gluon jets[8,16]. 

In addition to the above results which are consistent with our expectations on 
clan properties in parton showers, the model shows energy independent behavior 
(see Figure 9) by normalizing the average number of clans produced in a fixed 
rapidity interval \y\ < y c to the corresponding average number in full phase space, 
and by expressing this ratio as a function of the rescaled rapidity variable y* = 
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Vc/Vips- 



^W)= "%^ (79) 



This new regularity turns out to be stable for different choices of the parameter 
A. In Figure 9 a clean linear behavior is shown for the above ratio corresponding 
to the parameter value A=2. 

It should be pointed out that the asymptotic (log W — > oo) expression of the 
average number of clans for A = 2 corresponding to the interval log 2 < y c < 
logW — 2 log 2, which covers almost all the available phase space, is given by 

i v( !/c ^)^2 !/c i^ 2 + l0 «( l0 «^ l0 « 2) 



logW 

+ . n2 -^ 2 / c log2 + ^ c (log2) 2 -3 ? / c log21og(logW) + 

(log WO 2 L 3^3^ (8Q) 

+ y c log 2 log(log 2) - - (y c - log 2) 2 log(y c - log 2) + 



+ ^(y c + log2) 2 log(y c + log2) 



+ o 



1 



logW 



Eq. (80) can be obtained from eq. (A8) in the Appendix. To first order this 
formula shows that asymptotically the average number of clans in a fixed rapidity 
interval grows linearly with y c and is slowly decreasing with the energy of the initial 
parton. By including second order contribution, N(y c , W) continues to be linear 
in y c and to decrease with the energy but the slope of the curve is different from 
that at first order; this difference is quite remarkable (a factor 1.5) at low energies 
(below 100 GeV), but decreases at high energies. Notice that eq. (80) describes 
perfectly well the result for the average number of clans obtained by using the 
complete formula discussed in the Appendix in the same rapidity domain. 
The asymptotic expression (logVF — > oo) of n*(y*, W) is given by 



i+2io g (^; 

(1 

l + 21og(^f) 



+ M) 2 - 1] Mi - vt) - (i + vt? Mi + vt) , (81) 



. i//;:iog2 + (/y;:) 2 | 



1 + 2 log (YogT) log W log log W 
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Eq. (81) describes into 2% the result obtained by using the complete formula after 
inserting eq. (A8) into eq. (79). It should be pointed out that this behavior is 
maintained even at higher energies, indicating a scaling behavior in a wide energy 
range (50 500000 GeV). In the limit W -> oo, eq. (81) becomes 

^,W)^ + 0( 1 ^ W ) (82) 

It should be noticed that the slope of this curve is 10% less than that predicted 
by eq. (81), showing that the asymptotic regime is reached very late. 



Conclusions 

We have discussed a simplified model, based on essentials of QCD and lo- 
cal weakening of conservation laws for a shower (jet) originated by a parton of 
given maximum allowed virtuality and given rapidity. It generalizes previous re- 
sults to the case in which rapidity is added to virtuality evolution by assuming 
that clans are intermediate independent parton sources; clan production in this 
framework can be described as a Markoffian process. The new model, which we 
called Generalized Simplified Parton Shower model (GSPS), allows to determine 
clan multiplicity distribution in full phase space, which turns out to be a shifted 
Poissonian to be compared with the Poissonian behavior assumed in the standard 
interpretation of NB regularity. For clan multiplicity distributions in symmetric 
rapidity intervals, we obtain a combination of two Poissonian distributions. It is re- 
markable that the model predicts clan properties which show the same qualitative 
general trend observed both in experimental data and in Monte Carlo simulations. 
In addition it has been shown that the ratio N(y* yt ps , W)/N(W) is linear and 
energy independent when plotted as a function of rescaled rapidity variable y*. 

In this paper attention has been payed to the average number of clans pro- 
duced in a given rapidity interval by an initial parton of maximum allowed virtu- 
ality W, N(y c , W); the next step in our program is to calculate the corresponding 
average number of partons per clan in the same rapidity intervals, n c (y c , W). We 
decided to postpone this calculation in the framework of the GSPS to a forthcom- 
ing paper. 



Appendix 

In this Appendix we present a short discussion of the calculation of the in- 
tegrals given in eqs. (73) and (77) for A = 2. Although the integration of the 
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densities r^yolQy) an d P-^iUolW) might be considered straightforward, it is in- 
deed complicated by the appearance of very many terms. They correspond in case 
of r^(y \Qy) (eqs. (68) and (71)) to four regions of phase space which, according 
to Figure 10, we label with Roman letters. They are 

[a] y-\ogQ < y < y - log Q + log 2 

[b] y-log<2 + log2 < y < y 

(Al) 

[c] y < y < y + log<2-log2 

[d] y + logQ-log2 < y < y + logQ 



The densities r y (yo\Qy) are different in the regions [a],[b], [c],[d], and given by: 



r [a] 



(log 2W - y ) log Q + (log Qf log log Q+ (A2) 



r T \b] = 



(logW) 2 (log 2) 2 

2(logQ) 2 - (logQ) 2 log(log2W - yo - logQ) 

21ogQ 2(logQ) 2 log(logQ-log2) 



(\ogwy 



log 2 



(log 2)5 



2 log Q 2 log log Q 

H n — ^2 + log(2 log Q + y -\og2W)+ 

(log2) 2 

- log(log 2 W - y - log Q) + log(log Q - log 2) - log log 2 



rJc] = 



(\ogwy 



21ogQ 2(logQ) 2 log(logQ-log2) 
log2 



(log 2)5 



r T \d] 



2(logQ) 2 loglogQ 
+ (bi2p + lo g( lo S 2H/ - ^°) + 

- log(y - log 2W + log Q) + log(log Q - log 2) - log log 2 
2 2 



(log WO 2 (log 2) 2 



(y -\og2W) \ogQ+ 



- (log Q) 2 log(y - log 2W + log Q) + (log Q) 2 log log Q 



(A3) 



(A4) 



(A5) 



where the approximation defined in eq. (78) for y has been used. 

Accordingly, the density Pj,(yo\W) of eq. (73) results to be different in the 
above mentioned regions. Its integration over y variable leads to divide the in- 
tegration domain into five regions which we label by Greek letters in Figure 10, 
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i.e.: 

(a) -logW + log2 < y < -logW + 21og2 

(0) -logW + 21og2 < y < log2 

(7) log 2 < y < log W - log 2 (A6) 

(5) log W - log 2 < y < log W 

(e) log VP < y < log VP + log 2 

The analytic solution of the integral in eq. (77) turns out to be different for y c in 
the following different domains: 

Vc < log 2 
log 2 < y c < log W - 2 log 2 
log W - 2 log 2 < y c < log W - log 2 ( A7) 

log W - log 2 < y c < log W 

log W < y c < log W + log 2 

In order to give the flavor of the cumbersome and long calculations which lead to 
the final form of N(y c , W), whose asymptotic expression has been given in eq. (80), 
we show the result for y c in the domain log 2 < y c < log W — 2 log 2, i.e.: 
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- 2 / logTyiog2 (log2) 2 5y c ,log^ y c (logW) ; 
iV(y c , VP) = <^ + + 



(\ogW) 2 { 2 2 3 3 log2 

7y c log2 rt , rt 4y c (log2) 2 y 2 y 3 

\ + 2y c log W log2 + yc \* ' + f + 
6 3 2 3 log 2 

^ 2 , 4 y c (log H Q 3 
3 (log 2) 



4 y c (log HQ 3 2y c log 2 ^ 
2 y c log W 3(i og 2)2 g J log(log W - log 2) + 

(- log W + log 2) y c log(log 2) + 

^ (-( lo S W^) 2 - 4 log iy log 2 - 3 (log 2) 2 + 2 y c log 
d (log 2) 

4 y c log 2 - y 2 ) log(log W - y c )- 

(- log W + log 2 + y c ) 2 log(log W - log2 - y c ) , 



log ( log^ + log2-, c ) o§ + 8 o§ 2 + 3 o§ 2)2+ 

12 (log 2J Z 

2 log W (log 2) 3 + 5 (log 2) 4 - 8 y c (log W) 3 - 24 y c (log W) 2 log 2- 
18 y c log W (log 2) 2 - 14 y c (log 2) 3 + 12 (log W) 2 y 2 + 24 y 2 c log W log 2+ 
15 y 2 (log 2) 2 - 8 y 3 log W - 8 y 3 log 2 + 2 y c 4 ] + 

lQg g° g o ^ 2) t yc) ((log W) 2 + 4 log W log 2 + 3 (log 2) 2 + 2 y c log 

/i i o i 2\ n w . ^2 (-log2 + y c ) 2 log(- log 2 + y c ) 
4y c log2 + y c ) (logW + y c ) - + 

(log2 + y c ) 2 log(log2 + y c ) 



(log Ti" + log 2 + y c ) 4 log(log + log 2 + y c ) 
6(log2) 2 



(A8) 
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Figure Captions 



Fig. 1. The structure of the Simplified Parton Shower model (SPS) for parton 
production in virtuality Q and rapidity y. The degrading from the maximum 
allowed virtuality W is given by the product of eqs. (7) and (15). The splitting at 
(Q, y) is controlled by the product of eqs. (4), (10) and (14). Each dot represents a 
parton which further generates. Notice that here Qo + Qi < Q, \yo — y| < log 
Fig. 2. The structure of the Generalized Simplified Parton Shower model (GSPS) 
for clan production in virtuality Q and rapidity y. The production process is 
decoupled at each splitting both in virtuality and in rapidity (see eqs. (18) and 
(20)). Each blob represents a clan. Notice that here Qo < Q, Qi < Q, \yo — y\ < 
log^, \yi-y\ < log£. 

Fig. 3. Diagrams for the probability that an ancestor, splitting at virtuality 
Q and rapidity y gets virtuality Qo and rapidity yo after emitting N — 1 clans, 
^N(Qoyo\Qy) (solid line); dashed diagrams correspond to the interchange of the 
final ancestor and the last produced clan, i.e., to the probability to produce the 
(N — l) th clan with virtuality Qq and rapidity yo from an ancestor splitting at 
virtuality Q and rapidity y, which can be described by the same rN(Qoyo\Qy)- 
Fig. 4. Diagrams for the probability that an ancestor of maximum allowed virtu- 
ality W gets virtuality Qo an d rapidity yo after emitting N — 1 clans, Pn(QoVo\W) 
(solid line); dashed diagrams correspond to the interchange of the final ancestor 
and the last produced clan, i.e., to the probability to produce the (N — l) th clan 
with virtuality Qo and rapidity yo from an ancestor of maximum allowed virtuality 
W, which can be described by the same Pn(QoVo\W). 

Fig. 5. Phase space domain in virtuality and rapidity for an ancestor parton of 
maximum allowed virtuality W. Dashed lines represent the approximation used 
for explicit calculations (eq. (78)). The thick solid line in the lowest right corner 
corresponds to the contribution of one-parton shower, i.e., to the first diagram 
in Figure 4. In the central part of the figure we show the line corresponding to 
eq. (15). 

Fig. 6. Clan density p^(yo\W) for one shower at W = 50 GeV (dotted line), 100 
GeV (dashed line) and 500 GeV (solid line). 

Fig. 7. Clan density p^(yo\W) resulting by the addition of two back-to-back 
showers at W = 50 GeV (dotted line), 100 GeV (dashed line) and 500 GeV (solid 
line). 
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Fig. 8. Average number of clans N(y c: W) for a shower as a function of rapidity 
width y c at W = 50 GeV (dotted line), 100 GeV (dashed line) and 500 GeV (solid 
line). The symbols represent the values of the average number of clans in full 
phase space. 

Fig. 9. Normalized average number of clan n*(y*, W) for a shower as a function 
of rescaled rapidity y* = y c /yf ps at W = 50 GeV (dotted line), 100 GeV (dashed 
line) and 500 GeV (solid line). 

Fig. 10. Phase space domain in virtuality and rapidity for an ancestor parton of 
maximum allowed virtuality W. Different domains defined in eq. (Al) are labelled 
by Roman letters; different integration domains of eq. (A6) are labelled by Greek 
letters. 
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